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Probing the order-to-chaos region 
in superdeformed 151 Tb and 196 Pb nuclei with continuum 7-transitions 

S.Leoni 1 , G.Benzoni 1 , N.Blasi 1 , A.Bracco 1 , S.Brambilla 1 , F. Camera 1 , A.Corsi 1 , F.C.L.Crespi 1 , 
P. Mason 1 ' 2 , B. Million 1 , D.Montanari 1 , M.Pignanelli 1 , E.Vigezzi 1 , O.Wieland 1 , M. Matsuo 3 , 
Y. R. Shimizu 4 , P.Bednarczyk 5 , M. Castoldi 6 , D.Curien 7 , G. G.Duchene 7 , B. Herskind 8 , M. 
Kmiecik 5 , A. Maj 5 , W. Meczynski 5 , J.Robin 7 , J. Styczen 5 , M. Zieblinski 5 , K. Zuber 5 , A.Zucchiatti' 

1 Dipartimento di Fisica and INFN, Sezione di Milano, Milano, Italy 
2 Laboratori Nazionali di Legnaro, Padova, Italy 
3 Graduate School of Science and Technology, Niigata University, Japan 
4 Department of Physics, Kyushu University, Fukuoka, Japan 
5 The Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Krakow, Poland 

6 INFN sezione di Genova, Genova, Italy 
7 Institut de Recherches Subatomiques, Strasbourg, France and 
8 The Niels Bohr Institute, Copenhagen, Denmark 
(Dated: June 30, 2008) 

The 7— decay associated with the warm rotation of the superdeformed (SD) nuclei 151 Tb and 
196 Pb has been measured with the EUROBALL IV array. Several independent quantities provide 
a stringent test of the population and decay dynamics in the SD well. A Monte Carlo simulation 
of the 7-decay based on microscopic calculations gives remarkable agreement with the data only 
assuming a large enhancement of the B(E1) strength at low excitation energy, which may be related 
to the evidence for octupole vibrations in both mass regions. 
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The thermally excited, rapidly rotating atomic nucleus 
is an ideal laboratory for studying structure properties 
beyond mean field. In particular, the 7-decay which 
cools the hot compound nucleus at high angular momen- 
tum is a powerful tool for investigating the transition 
from the chaotic regime around the particle binding en- 
ergy to the cold, ordered system close to the yrast line 
[Hi M El- This topic, mostly investigated in normally 
deformed (ND) nuclei [l|, [H, |j, [f| 0, @], becomes partic- 
ularly interesting and challenging, both experimentally 
and theoretically, in the case of very elongated systems. 
The study of superdeformed (SD) nuclei requires, in fact, 
high selectivity and high statistics experiments, in order 
to focus on the small fraction of the 7-decay (« few %) 
which finally gets trapped into the SD well and populates 
discrete rotational bands. At present, there is only par- 
tial experimental information available for few SD nuclei 
d, [TO, EH E3] > whose interpretation is based on models 
which arc cither rather schematic [9|, [ll|, ll2j or depend 
on several parameters [1, [l(| • 

The aim of this letter is to make a substantial step 
forward in the understanding of the population and de- 
cay of warm SD nuclei in the mass regions A=150 and 
A=190, where superdeformation is a well established 
phenomenon, by a comparative and comprehensive study 
of the nuclei 151 Tb and 196 Pb. For the first time, sev- 
eral independent observables are extracted from quasi- 
continuum 7-coincidence spectra. This provides strong 
experimental constraints on the dynamics of the 7-decay 
flow and of the tunneling through the potential barri- 
ers between SD and ND excited states, over a wide spin 
range. A key point of this work is the development of a 
new Monte Carlo simulation of the 7-cascades in the SD 



well. This makes use of discrete levels calculated micro- 
scopically with the cranked shell model. In addition, the 
coupling between excited SD and ND configurations is in- 
cluded, based on penetration probabilities obtained from 
a microscopical model [13]. It is anticipated that such 
a detailed study of the order-to-chaos region reveals the 
influence of nuclear structure effects. In particular, an 
enhancement in the El strength around 1-2 MeV, which 
could be related to octupole vibrations, is found to play 
an important role. 

The experiments were performed at the Vivitron in 
Strasbourg (F) with the EUROBALL IV array [H|. The 
reaction used were 130 Te( 27 Al,6n) 151 Tb (at 155 MeV) 
and 170 Er( 30 Si,4n) 196 Pb (at 148 MeV). In both cases, 
a stack of two self-supporting targets (with a total thick- 
ness of ~ 1 and 1.2 mg/cm 2 , respectively) were employed. 
In the Tb case, the full Ge ball was used, while in the 
Pb experiment the low efficiency Ge detectors in the for- 
ward hemisphere were replaced by 8 large volume BaF2 
scintillators to measure high energy 7-rays from the giant 
dipole resonance (GDR). In both cases an InnerBall of 
BGO detectors was used to determine the 7-multiplicity 
of every event. The SD nuclei 151 Tb and 196 Pb were 
investigated in detail by discrete 7 spectroscopy [l5j . 
In both cases discrete transitions linking the SD yrast 
band to ND configurations were identified, and in the 
196 Pb case it was possible to determine the spin and ex- 
citation energy of the SD yrast band [l(| [13] • In this 
work we focus instead on the unresolved excited SD ro- 
tational bands, which form ridge and valley structures 
in 7-coincidence spectra. The two data sets have been 
sorted into a number of 7 — 7 matrices in coincidence 
with low-lying ND transitions of the isotope of interest 
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( 151 Tb or 196 Pb), named Total, and in coincidence with 
the SD yrast band of each nucleus (named SD-gated). 
In addition, rotational planes (ROT-Plane) , namely ma- 
trices obtained selecting triple coincidences with the re- 
quirement x + y = 2z, x, y and z being the energies of the 
7-rays, were also constructed in coincidence with the nu- 
cleus of interest [l8j , in order to enhance the sensitivity 
to rotational correlation in the 7-cascades. In all cases 
a condition on high-fold events (F > 25 for 151 Tb and F 
> 10 for 196 Pb) was imposed, to better focus on high- 
multiplicity cascades. For each matrix the background 
was subtracted by scaling the ungated spectrum by a 
factor corresponding to the P/T of the ga ting transitions 
and further applying the COR procedure [l9j ] in order to 
eliminate the remaining uncorrelated background. In all 
spectra cuts perpendicular to the main diagonal reveal 
the existence of ridge structures extending over a broad 
interval of transition energies (corresponding to the spin 
regions 40 to 60 % and 10 to 40 h for 151 Tb and 196 Pb), 
with a spacing between the two most inner ridges equal to 
2AE 7 = 8h 2 /5;W ~ 100 and 80 keV respectively, where 
3( 2 ) is the moment of inertia of the SD yrast band in 
each nucleus. Figure [T] shows examples of such cuts for 
the wide spin intervals w 48-56 Ti ( 151 Tb) and w 20-34 
% ( 196 Pb). In the case of the Total and SD-gated spec- 
tra (a), b) and e), f)) the dotted lines represent the full 
intensity distributions, while the solid lines show the in- 
tensity remaining after the subtraction of all known ND 
and SD discrete peaks, making use of the RADWARE 
package [20|]. The ROT-planes spectra (c) and d)) are 
sum of cuts taken in between the SD yrast peaks. In this 
case, the innermost ridge corresponds to the second ridge 
in the Total and SD-gated matrices, and the associated 
spacing is equal to 4 Ai? 7 . 

In first place, the intensity of the ridge structures has 
been evaluated and compared to the SD yrast, as shown 
by symbols in figure [5] In both nuclei the total intensity 
of the 1 st and 2 nd ridge (panels a), b) and c), d)) is up to 3 
times larger than the SD yrast population at the plateau 
(which collects 2% and 1.3% of the total decay-flux, re- 
spectively), pointing to the existence of several discrete 
unresolved SD bands decaying directly into the low defor- 
mation minimum. The more selective analysis of the 1 st 
ridge in coincidence with the SD yrast (panels e) and f)) 
shows, in fact, that only a fraction of the total population 
of the discrete excited bands is finally collected into the 
SD well, corresponding at most to 40% and 80% of the SD 
yrast intensity in 151 Tb and 196 Pb, respectively. Finally, 
the symbols in panels g) and h) give the total intensity 
of the rotational decay in coincidence with the SD yrast, 
including also the contribution from fragmented/damped 
bands at higher excitation energies. This is obtained by 
the analysis of the E2 bump observed in one dimensional 
spectra, in coincidence with the SD band, which in the 
case of Pb show the presence of contaminants (probably 



Mi's) in the spin region 20-30 h [21(. It is found that 
the intensity of the SD-gated ridges accounts for almost 
the entire E2 strength up to spin 55 fi in 151 Tb and 35 
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FIG. 1: Projections perpendicular to the main diagonal of 7- 
coincidence matrices of 151 Tb (left) and 196 Pb (right), for the 
energy interval 1150-1280 keV (i.e. w 48-56 h) and 430-680 
keV (i.e. 20-34 h), respectively. Panels a), b) show cuts on 
the Total matrices, c), d) from the rotational planes and e), 
f) from the SD-gated spectra (see text for details). 



h in 196 Pb, corresponding to the plateau regions of the 
SD yrast, while at higher spins the ridge intensity drops 
to less than half of the E2 bump, indicating that most 
of the feeding of the SD yrast comes from fragmented 
rotational bands in the region where rotational damping 
largely dominates [l[ . This differs from the peculiar case 
of 194 Hg, were an exceptionally narrow ridge, exhausting 
nearly all E2 decay strength has been interpreted as a 
signature of ergodicity in a nuclear system [l2j . 

Additional information on the dynamics of the 7-decay 
in the SD well and on the coupling with the ND states 
can be obtained estimating the number of discrete bands 
populating the ridge structures from the analysis of the 
fluctuations of the counts in the 7-coincidence matrices 
[J [|J . The results are shown by symbols in figure for 
151 Tb (left) and 196 Pb (right). In both cases a rather 
large number of SD discrete bands (up to more than 30) 
are found to populate the Total ridges (a) and b)), half 
of which extend at least over 3 consecutive transitions, 
as follows by the analysis of the second ridge (c) and d)). 
In addition, the study of the SD-gated ridges (e) and f)) 
show that no more than half of the discrete excited bands 
arc feeding into the SD yrast, supporting the analysis of 
the ridge intensities, previously discussed. 

The experimental findings have been compared with 
predictions based on a Monte Carlo simulation of the 7- 
decay flow of 151 Tb and 196 Pb. The code is an extended 
version of MONTESTELLA [H and simulates the 7- 
decay of the excited rotating nucleus from the residual 
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FIG. 2: Population of the excited rotational bands in the 
SD well of 151 Tb (left) and 196 Pb (right), relative to the cor- 
responding SD yrast (normalized to 100% at the intensity 
plateau and shown by circles in the bottom panels). The 
symbols show the experimental data: the intensities of the 
I s * and 2 nd ridge are given in panels a), b) and c), d); panels 
e), f) refer to the SD-gated ridges, while panels g), h) show 
the total E2 population in coincidence with the SD band, as 
extracted from the E2 bump (see text). Predictions from sim- 
ulations of the 7-decay assuming a standard (enhanced) El 
strength (see top panels i) and 1)) are shown by dotted (solid) 
lines. 



entry distribution towards the SD or the ND yrast line. 
It is based on the competition between E2 collective and 
El statistical transitions in both minima and on a tun- 
neling probability across the barrier separating the two 
potential wells, accordingly to the equations of ref. [23| . 
In the SD well energy levels, E2 transition probabilities 
and potential barriers are microscopically calculated ac- 
cording to the cranked shell model of ref. [HI, |24|, in- 
cluding, in the case of Pb, the phenomenological adjust- 
ment of the mass parameter C mass = 2.7 to reproduce 
the decay-out spin of the SD yrast. At excitation en- 
ergies higher than the ones covered by the microscopic 
states (i.e. U> 4.5 and 3.5 MeV in 151 Tb and 196 Pb, re- 
spectively) the rotational decay in the SD well proceeds 
through a continuum of states, governed by level densi- 
ties, E2 strengths (with a width T rot ~ 100, 70 keV, re- 
spectively) and barriers extrapolated from the region of 
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FIG. 3: The number of discrete excited SD bands obtained 
from the analysis of the count fluctuations of the ridge struc- 
tures of 151 Tb (left) and 196 Pb (right). Panels a), b) and 
c),d) refer to the I s * and 2 nd ridge, e), f) to the SD-gated 
ridge. Symbols refer to data, solid lines to the Monte Carlo 
simulation. See the text for the dashed and dotted lines. 



microscopic levels. The 7-decay in the ND well is instead 
described schematically, as done in the simulations of ref. 
B lipt. with a density of states calculated accordingly to 
ref. [25l| . Once the entry distribution of the 7-decay is de- 
fined, the use in the SD well of microscopically calculated 
quantities leaves no room for adjustable parameters, ex- 
cept for the El transition probability, which is described 
in terms of the Lorentzian strength function of the gi- 
ant dipole resonance, with an hindrance factor tuned to 
reproduce the intensity of the SD yrast (as in the code 
MONTESTELLA [H, The residual entry distri- 

butions for the 7-decay of 151 Tb and 196 Pb have been 
calculated in two steps. Firstly, we have used the heavy 
ions-collision model of ref. [26[ , which provides the fusion 
cross section of the compound nuclei 157 Tb and 200 Pb. 
Secondly, we have determined the neutron evaporation 
by the Monte Carlo version of the code CASCADE [27[, 
obtaining angular momentum and excitation energy dis- 
tributions of the 7-decay with a maximum at I = 60 and 
40 h and U = 8.3 and 7.7 MeV above yrast, for 151 Tb 
and 196 Pb, respectively. 

The Monte Carlo simulation provides predictions for 
all the observables shown in figure and [31 including, 
for the first time, the intensity and number of discrete 
bands in coincidence with the SD yrast, which depend 
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strongly on the competition between E2 and El transi- 
tions along the cascades. For both nuclei the model re- 
produces rather well the more inclusive quantities, such 
as the SD yrast intensity and the population of the dis- 
crete excited bands forming the I s * and 2 nd ridge (dotted 
lines in figures 2 a), b) and c), d)). On the other hand, 
it largely underestimates the E2 intensity observed both 
in the ridge and E2 bump in coincidence with the SD 
yrast (dotted lines in panels e),f) and g),h)). These ob- 
servables, although referring to the warm rotation in the 
transition region between order and chaos, are in fact 
very sensitive to the delicate balance between E2 and El 
transitions at low excitation energy, where nuclear struc- 
ture effects play a major role. It is worth noticing that 
in both A=150 and A=190 mass regions, experimental 
evidence has been found for octupole vibrations in su- 
perdeformed nuclei ( 152 Dy [28], 196 Pb [H, 190 Hg [H), 
resulting in strongly enhanced El transitions, linking the 
excited SD bands to yrast. The B(E1) values are in fact 
of the order of 10~ 4 W.u., namely 1 to 2 orders of mag- 
nitude larger than expected from the tail of the GDR, 
in agreement with various theoretical models [3ll |32| . 
We have therefore performed simulation calculations for 
151 Tb and 196 Pb making use of a GDR strength with an 
extra component giving rise to an enhancement in the 
B(E1) values between 10 to 100, in the excitation energy 
region around 1.5 and 1 MeV, respectively, i.e. consistent 
with the energies of the El enhanced transitions experi- 
mentally observed (see figure 2 i) and 1)). The obtained 
results (solid lines in figure 2) are now much better in 
agreement with the data, especially for the E2 strength 
in coincidence with the SD yrast (panels e), f) and g), 
h)). The number of paths N pat h — l/Y)wf, where u>i in- 
dicates the relative population wt of each path, is also in 
good agreement with the data (solid lines in Fig. [3]). 
This quantity represents a crucial test for various as- 
pects of our model. N pat h gives the number of discrete 
bands which are populated sufficiently strongly and it is 
determined by the level density, by the onset of damp- 



ing (which takes place at U ~ 2 MeV in 196 Pb and at 
U ~ 3 MeV in 151 Tb) and by the population mechanism. 
The latter favours low- lying bands around 1-1.5 MeV, 
not only because they tunnel less easily through the po- 
tential barrier but also because they receive stronger El 
feeding. We note that N pat h is much lower than the num- 
ber of discrete bands Nb an d obtained from the cranked 
shell model ignoring the flow (dotted lines in Fig. tak- 
ing into account the depopulation of the bands due to the 
tunneling (as estimated in ref. [l3j]) decreases Nb an d at 
low spins but still largely overestimates the data (dashed 
lines). This is different from the case of ND nuclei, where 
discrete bands are more equally populated and N pat h is 
closer to N band ^2J. 

In summary, we have presented a study of the warm 
rotation in the SD nuclei 151 Tb and 196 Pb, providing 
a series of independent quantities which test our under- 
standing of the 7-decay in the SD well over the whole spin 
range, i.e from the feeding region, down to the decay-out 
into the ND minimum. A newly developed Monte Carlo 
code, based on microscopic quantities, is able to repro- 
duce all the observables except the E2 strength gated on 
the SD yrast line. Agreement with this quantity can be 
obtained only increasing the El strength at low excita- 
tion by almost two orders of magnitude, as compared to 
the standard Lorentzian parametrization for SD nuclei. 
This may be related to the observation of strong El tran- 
sitions between discrete bands, observed in both regions 
and associated with octupole vibrations. Therefore, the 
study of the warm rotation in the transition region be- 
tween order and chaos not only is instrumental in shed- 
ding light on nuclear structure properties beyond mean 
field, but it can also be considered as a tool to unravel 
microscopic features of the 7-decay. 
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